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Abstract 

We compare the transverse momentum (pr) distribution of inclusive light-charged- 
particle production measured by the CDF Collaboration at the Fermilab Tevatron 
with the theoretical prediction evaluated at next-to-leading order in quantum chro- 
modynamics (QCD) using fragmentation functions recently determined through a 
global data fit. While, in the lower px range, the data agree with the prediction 
within the theoretical error or slightly undershoot it, they significantly exceed it 
in the upper pr range, by several orders of magnitude at the largest values of px, 
where perturbation theory should be most reliable. This disagreement is too large 
to be remedied by introducing additional produced particles into the calculation, 
and potentially challenges the validity of the factorization theorem on which the 
parton model of QCD relies. Clearly, a breakdown of the factorization theorem, 
being a fundamental property of QCD, would be extremely difficult to understand. 
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Inclusive single-hadron production processes in e + e~ , ep, and pp(p) reactions can be 
well described [U [2] by cross section calculations that depend on universal fragmentation 
functions (FFs), consistent with the validity of the factorization theorem at CERN LEP, 
BNL RHIC, and DESY HERA energies (for a recent review, see Ref. [3]). However, in 
order to prepare for incoming data at much higher energy from the CERN LHC, from the 
late 2009 run and/or from future runs, this validity needs to be tested also at the Tevatron. 
Such tests were performed in Ref. p] , using CDF data for which p T reached 10 GeV jl] and 
also UAl and UA2 data from CERN SppS runs at sub-TeV energies. Very recently, new 
data have been obtained by the CDF Collaboration [5] for which p? reaches 150 GeV. Such 
data provide a unique opportunity to test the factorization theorem applied to the final 
state in a high-p^ regime of inclusive hadron production that has never been accessible 
before, where we expect perturbation theory in the QCD-improved parton model to work 
perfectly, and these data are the focus of this Letter. We note that agreement between 
theory and recent CDF measurements of inclusive jet production [B] has been obtained 
at much higher px values, which confirms the inital-state factorization theorem for these 
data. Such measurements originate from the same physical processes as the lower-p^ 
hadron production measurements do, since the produced hadron carries just a fraction of 
the jet's energy. 

We present here the first rigorous interpretation of the CDF data of Ref. [5] in the 
context of perturbative QCD. The experimental analysis resorted to the Monte Carlo 
event generator PYTHIA, which failed to give useful results for p? > 50 GeV. We use 
the latest FF sets [9j EJ [7J extracted from global data fits that included accurate primary- 
quark-tagged measurements from e + e~ reactions. Unfortunately, e + e~ reaction data do 
not sufficiently constrain the gluon FF, because it only enters the calculation at next- 
to-leading order (NLO), while it enters the calculation of pp reactions at the same order 
as the quark FFs. The two most recent FF sets, AKK08 [7J and DSS [5], also used pp 
reaction data from RHIC, but the gluon FFs in both cases differ substantially, suggesting 
that the RHIC data used in those fits did not significantly improve the constraints on the 
gluon FFs. The FF errors in both sets were not determined. Therefore, the best estimate 
of the error on the predictions for CDF data at present is given by the general spread in 
the calculation using different sets. Besides the AKK08 and DSS FF sets, we also include 
the HKNS one [9], which was obtained from a fit only to e + e~ data. We assume that 
the unidentified charged particle in these CDF measurements is a light charged hadron. 
That is, we take the sum of the n^^, K ± and p/p production cross sections, and assume 
the possible contamination with other charged particles to be negligible. 

Our calculations are performed to NLO in the modified minimal-subtraction (MS) 
renormalization and factorization schemes. We set the renormalization scale to be \x = kpr 
and the factorization scale to be /i/ = kjpT, where k — kf — 1 unless otherwise stated. For 
the initial protons, we use the parton density function (PDF) set CTEQ6.6M |10j . Since 
the produced-hadron mass was accounted for only in the AKK08 analysis, we account for 
it in the same way here for calculations using the AKK08 FF sets, but set it to zero when 
using the DSS and HKNS sets. In any case, neglecting the hadron masses makes very 
little difference, as we shall see later. We refer to the above calculation as the AKK08 
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default one (labeled "AKK08(default)" in our figures). 

The comparison with the recent data [5] collected in Run II (with center-of-mass energy 
y/s = 1.96 TeV) in the central region of the CDF detector (with pseudorapidity \r)\ < 1) 
and older CDF data jl] from Run I (y/s = 1.8 TeV, < 1) is shown in Fig. [TJ The 
most striking observation to be made is the disagreement at px > 50 GeV (see left plot). 
Looking closely at the experimental line shape, there appears to be a transition to a softer 
slope, corresponding to a fall-off with a lower power in 1/px, at about pt = 20 GeV, which 
eventually leads to a systematic departure from the theoretical prediction. In fact, fitting 
the function A/p T to the CDF data yields A = 71.50 and n = 6.52 with xV d -°- f - = °- 886 
for 5 GeV < p T < 20 GeV and A = 0.030 and n = 3.97 with x 2 /d.o.f. = 2.18 for 
20 GeV < px < 150 GeV. We evaluate the theoretical error, indicated by the shaded 
region labeled "AKK08(Min/Max)" in Fig. [1] and later plots, by plotting, for each px 
value, the maximum and minimum cross sections in the range 1/2 < k, kf < 2. To 
achieve this, we use a 17x17 grid of points equidistant in ln/c(/) ranging from — In 2 to 
In 2. As Fig. [TJ shows, the theoretical error calculated in this way is not nearly sufficient 
to accommodate the recent CDF data at large px- Within this error, the calculation in 
the region px ^ 40 GeV is consistent with most of the data. In the region px ^ 20 GeV, 
the data lie in the lower range of the theoretical error, and in fact favour a large value of 
k = kf, as can be seen by the dashed curve in Fig. [TJ which corresponds to k = kf = 4. 
The error is largest at small px = 0(1), where both a s (/i) and the partonic cross sections 
become very large, and for this reason we limit our discussion to those data for which 
Px > 2 GeV. 

In Fig. [2j the same calculation is performed with the DSS [11] and HKNS [9] FF sets. 
We note that FFs for unidentified charged hadrons were also obtained by the authors of 
Ref. [11] . However, the difference between the calculation using those FFs and the one 
using the sum of the DSS FF sets for tt^, K^, and p/p production (which respectively make 
up about 75, 15, and 10% of the sample) turns out to be negligible over the whole px range 
considered. Assuming that the difference between the calculations using the AKK08, DSS, 
and HKNS FF sets gives an estimate of the error due to the uncertainties in the FFs, we 
conclude that the latter are not nearly sufficient to account for any discrepancy between 
theory and the new CDF data [5] . We stress that all these FF sets are tightly constrained 
in the very region of the (z,/if) parameter space that is probed by our predictions, in 
particular where the discrepancy occurs, so that no extrapolation error arises. In fact, 
the average value of z lies in the range 0.44 ^ (z) <^ 0.65 for 2 GeV < \if = px < 150 GeV. 
The wealth of data from TRISTAN, LEP1, SLC, and LEP2, which populate the range 
58 GeV < [if = y/s < 191 GeV, are quite confining in this intermediate z range. We 
conclude that the error due to the uncertainties in the FFs is small against the overall 
one. 

In Fig. [3j we investigate the effect of the produced hadron's mass mh, by comparing 
our default prediction with the one obtained by setting all particles' masses to zero. The 
difference also gives an estimate for the size of any neglected small-p^ effect of a similar 
order of magnitude, namely 0(rn\/p T ), such as higher twist. Clearly, such effects are 
negligible relative to the theoretical uncertainty shown in Fig. [TJ 
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As mentioned above, the gluon FF is somewhat less well constrained than the quark 
FFs. In order to determine whether the discrepancy with the recent CDF data [5] could 
be due to large errors on the gluon FF, we also show in Fig. [3] the contribution coming 
from quark fragmentation only (labeled "no gluon"). We conclude that a reasonable 
modification to the gluon FF may improve the description of the low-j>r data, but certainly 
not at the largest px values where gluon fragmentation becomes negligible. 

The only other non-perturbative quantities in our analysis besides FFs are PDFs. Also 
these are tightly constrained by global analyses in the range of x and fif probed by the 
CDF data. In order to illustrate this, we study in Fig. H] the variation due to the choice 
of PDF set, by comparing our default result based on the CTEQ 6.6M PDF set \±0\ with 
the calculations using the MSTW2008 p2] and HERAPDF0.1 [13] ones. All these PDF 
sets have been constrained in very different ways. Specifically, the HERAPDF0.1 PDF 
set [13] was obtained from a combined fit to HERA data, which puts strong constraints 
on the gluon PDF at low x. Clearly this type of error is the least significant of all. 

Finally, we compare with the recent accurate data from STAR [T4] (y/s = 200 GeV, 
\y\ < 0.5) in Fig. [51 and find perfect agreement with our default prediction, well within the 
theoretical error. By contrast, the CDF data from Runs I |4j and II [5] tend to undershoot 
the theoretical predictions in the same px range (see Fig. [I]). A possible explanation for 
the worse agreement at the larger a/s values is that perturbation theory eventually fails 
if the fraction of available momentum z = 2px cosh y/\/s taken away by the produced 
hadron becomes too low. 

In conclusion, the pt distribution of inclusive charged-hadron hadroproduction re- 
cently measured by the CDF Collaboration [5] at the Tevatron was confronted here for 
the first time with rigorous NLO predictions relying on the factorization theorem of QCD, 
which, unlike MC event generators, have no left-over adjustable parameters. Surprisingly, 
a spectacular disagreement was discovered, by up to 3 orders of magnitude, at the largest 
Pt values, where perturbation theory is expected to work most reliably. In fact, the x 2 P er 
degree of freedom turned out to be as large as 7.3 for the last 4 data points. Perfoming a 
careful error analysis, we estimated the overall uncertainty at the largest pt values to be 
about ±30%, being mainly due to unknown corrections beyond NLO manifesting them- 
selves in residual scale variations. In fact, the theoretical uncertainty actually reduces with 
increasing value of pr simply because of asymptotic freedom causing the strong-coupling 
constant to fade, so that the high-p^ region is where we would expect the agreement with 
the data to be the best. Being at most of 0(Aqcd/pt), higher-twist effects are bound to 
be negligibly small there. The only large logarithms that need to be resummed are those in 
ct s (n), the PDFs, and the FFs, which are duly taken care of by the renormalization group 
and evolution equations, respectively. The partonic cross sections are devoid of large log- 
arithms requiring resummation, which is ensured by the scale choice jj, — {if — 0(px)- 
Because (z) ~ 0.5, no large logarithms of the types In z or ln(l — z) need to be resummed, 
and resummation may not even be appropriate. Because even the highest px bin is far 
away from the edge of phase space, no high-p^ resummation is required either. A different 
philosophy of error analysis may lead to a somewhat different result, but certainly not 
to an error of O(10 5 %), required to remove the disagreement. The latter is even far too 
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large to be remedied by introducing new heavy virtual particles with resonating propaga- 
tors or charged particles in the final state. Some more drastic modification either to the 
theoretical preduction or to the experimental data is required. We do not claim that the 
discrepancy discovered by us actually challenges the factorization theorem, nor that the 
data is wrong, nor that new physics is at work. We do not even wish to speculate about 
possible problems in the experimental analysis, such as the suppression of the cosmic-ray 
background, etc. We only state that there is a discrepancy, whatever its reason might be. 
We recall that no such discrepancy was observed in the inclusive single-jet production data 
of the same experiment [6] at considerably higher values. High-p^ data from the LHC 
is needed to verify whether this apparent final-state factorization breaking is a genuine 
phenomenon. We note that the recent inclusive charged-hadron production data from the 
ATLAS and CMS Collaborations [15] at the LHC are limited to the region px < 20 GeV. 
Unfortunately, they are presented as yields differential in p T and cannot be converted to 
cross section distributions because the conversion factors are not presently known. 

This work was supported in part by BMBF Grant No. 05H09GUE, by DFG Grant 
No. KN 365/5-3, and by HGF Grant No. HA 101. 

Note added: We were informed by the CDF Collaboration that they are currently 
reanalyzing their published data [5]. In the meantime, four more papers [16J on the same 
subject have appeared. 
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Figure 1: Top left: The CDF data of Ref. [5] for inclusive unidentified-charged-particle 
production is compared with NLO QCD predictions evaluated using the AKK08 FF sets 
[7]. The shaded region labeled "AKK08(Min/Max)" indicates, for each value, the 
maximum and minimum values of the cross section within the ranges 1/2 < k < 2 and 
1/2 < kf < 2. Also shown is the result for k — kf — 4. Top right: As in the top left 
plot, except for the older CDF data of Ref. jl]. Each lower plot is identical to the one 
above it, except that all results are divided by the central NLO predictions for clarity. 
(The latter hence correspond to unity.) 
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Figure 2: As in Fig. [T] (bottom left). Also shown is the calculation performed using the 
DSS HO and HKNS [S] FF sets. 
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Figure 3: As in Fig. [T] (bottom left) for the data. The AKK08 calculation is performed 
with all hadrons' masses set to zero (labeled "m^ = 0" ) and again without the contribution 
from the evolved gluon FF ( "no gluon" ) . 
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Figure 4: As in Fig. [TJ (bottom left) for the data. The AKK08 calculation is performed us- 
ing the PDF sets of Ref . [Lj (labeled "MSTW") and Ref. P] (labeled "HERAPDF0.1"). 
The data for which > 50 GeV lie far above the calculation. 
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Figure 5: The STAR data of Ref. [H] for inclusive charged-hadron production is compared 
with NLO QCD predictions evaluated using the AKK08 FF sets [7]. 
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